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ABSTRACT 



We study pair production of leptoquarks in e+e^ annihilation at linear 
colliders. A detailed simulation including beamstrahlung and initial state 
radiation, leptoquark decay and hadronization, as well as detector smearing, 
is performed. Discovery limits are estimated for center-of-mass energies of 
500 and 800 GeV. The prospects for determining masses and couplings of 
leptoquarks are also investigated. 
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1 Introduction 



Leptoquarks appear in extensions of the standard model involving unification, techni- 
color, compositeness, or i?-parity violating supersymmetry. Hence, the search for lepto- 
quarks is an important task at present and future high energy experiments. Furthermore, 
leptoquark pair production leads to novel experimental signatures involving fixed-mass 
lepton-jet systems, and thus provides an interesting example for detector studies. 

In e^e~ annihilation, pair production of scalar leptoquarks by s-channel 7 and Z 
exchange is uniquely determined by gauge symmetry if one assumes minimal couplings, 
while vector leptoquarks may also possess anomalous couplings. In addition to s-channel 
production, there are t-channel processes involving unknown Yukawa couplings. Simi- 
larly, pair production of leptoquarks at hadron colliders is determined by the QCD gauge 
couplings with contributions from anomalous and Yukawa couplings. In contrast, in 
lepton-nucleon scattering the dominant s-channel leptoquark production, as well as t- 
and w-channel scattering, proceeds via Yukawa couplings only. 

In the generally adopted framework described in Ref. [1], the Yukawa couplings are 
taken to be dimensionless and SU{3) x SU{2) x U{1) symmetric. Moreover, they are 
assumed to conserve lepton- and baryon-number in order to avoid rapid proton decay, to 
be non-zero only within one family in order to exclude FCNC processes beyond the CKM 
mixing, and chiral in order to avoid the very strong bounds from leptonic pion decays. 
The allowed states can be classified according to spin, weak isospin and fermion number. 
They are summarized in Table 1. 

The leptoquark masses and couplings are constrained by high-energy data. Direct 
searches for leptoquarks have been performed at the Tevatron, at HERA and at LEP. 
Recently, both experiments, CDF and DO, have improved their mass limits for scalar 
leptoquarks considerably. DO excludes first generation leptoquarks with masses below 
225 GeV assuming a branching ratio B^q = 1 for decays into electrons [2], whereas CDF 
quotes a limit of 213 GeV [3] (all mass limits are at 95% CL). For branching ratios less 
than one, the limits are weaker, e.g., M > 176 GeV for Beq — 0.5 [2]. The bounds on 
vector states are even stronger: 298 GeV for B^g — 1 and 270 GeV for B^g — 0.5 [4] . The 
corresponding bounds on second and third generation scalar leptoquarks are M > 184 
GeV for = 1 and M > 98 GeV for Brq = 1, respectively [5]. The LHC is expected 
to reach the TeV range. Mass limits obtained at HERA depend on the Yukawa couplings 
Xl,r- They range from 207 to 272 GeV for Xl,r — e [6] where e is the electromagnetic 
coupling strength . The precise mass limit also varies with the type of leptoquark specified 
in Table 1. The above limits are lowered by about 50 GeV if Xl,r = 0.1. Finally, the 
most stringent but Ai^R-dependent limit from LEP originates from the search for single- 
leptoquark production and excludes masses below 131 GeV assuming Xl^r > e [7]. The 
mass bounds from leptoquark pair production roughly reach \/i/2 [8], ^/s being the 
center-of-mass energy, and are thus much weaker than the Tevatron bounds. 

Indirect bounds on Yukawa couplings and masses can be derived from low-energy data 
[9]. For chiral couplings, the most restrictive limits come from atomic parity violation 
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State 






Channel 


Limits 


9i{s) 


(7tot(s)[fb] 
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Si 


-'/'So 


9l 

9r 
-9l 


(^rUr 


gL < 0.06 

gR < 0.1 


0.236 


6.15 


2 


Si 


-'/'So 
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e^dR 
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3.77 


98.2 






+V^Si 
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S3 


-4/3^^ 


-9l 
-9l 
-V2gL 
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gL < 0.09 
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6.05 


6.15 
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R2 


-2/3 Q 
' '-'1/2 

-'/'Si/2 


9l 
-9r 

9l 
9r 


VeUL 
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^R^R 


gL < 0.1 
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DO.O 
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-'/'Si/2 


Ot 

9l 
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gL < 0.1 


1.06 
1.51 


27.6 
39.5 
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— 1/3t ^ 
' Vi/2 

-^/'Vu2 


9l 
9r 
9l 
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l^edR 

^rUl 
eldR 

Cr(Il 


gL < 0.09 

Qr < 0.05 


1.56 

3.84 


365 

895 


7 


V2 


+2/3 f/ 
-'/'Vy2 


9l 
9l 


J^eUR 

eluR 


gL < 0.09 


i.ol 
1.06 


3o3 
247 


8 


Ui 


-^I^Vo 


Qt 
9r 
9l 


T^eUR 


gL < 0.05 
gR < 0.09 


0.942 


222 


9 


Ui 


-5/3^ 


9r 


CrUl 


gR < 0.09 


5.89 


1370 






+1/3^^ 


V2gL 


VedR 




4.02 


942 


10 


Us 


-5/3y^ 


-9l 
9l 
V2gL 


eldR 

T^eUR 

eluR 


gL < 0.04 


0.942 
7.67 


222 
1790 



Tabic 1: Properties of leptoquarks. Columns 2 and 3 show the notations of Refs. [1] and 
[12], respectively. The upper half of the Table refers to scalars, the lower to vectors. The 
lower indices denote the weak isospin, multiplicity 2T + 1 or total isospin T . The upper 
index is the electromagnetic charge. The Yukawa couplings gL,R are as defined in Ref. 
[1], while \l,r is used as a generic symbol. Column 5 specifies the production and decay 
channels, while the low-energy limits on the Yukawa couplings from Ref. [9] for M = 200 
GeV are given in column 6. The last two columns quantify the effective coupling gi(s) 
defined in Eq. (6) and the total cross section atot{s) given in Eqs. (11, 12) taking y/s = 500 
GeV, M — 200 GeV, gL — gR — 0, and including corrections due to beamstrahlung and 
initial state radiation. 
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and lepton and quark universality. The maximum allowed couplings for M — 200 GeV 
and first generation leptoquarks are given in Table 1 [10]. 

The present work extends previous studies [11] in several ways: 

(i) complete simulation of leptoquark production and decay 

e+e~ ^ ^ Iiqil2q2 or hqil2q2, (1) 
taking into account beamstrahlung, initial-state radiation and quark hadronization; 

(ii) simulation of the main background processes 

e+e" ^ WW, ZZ, (2) 

e+e~ hkqsQA, (3) 
e+e~ tt; (4) 

(iii) inclusion of acceptance and smearing efi^ects for different detector models. 

After the presentation of calculational details in section 2 and experimental considerations 
in section 3, the results on detection efficiencies and sensitivity limits are discussed in 
section 4. 



2 Theoretical Framework 



In e^e~ annihilation, leptoquarks can be pair-produced via photon and Z boson s-channel 
exchange and via i-channel exchange of quarks. The couplings of leptoquarks to gauge 
bosons are given in Ref. [13]^, the couphngs to fermions in Ref. [1]. With these couplings, 
the differential cross section for the production of scalar or vector leptoquarks of mass M 
is given by [13] 



da 



dcos9 



~8~ 



N { sin^ 9 



+ 5o + 



t 



(5) 



where t = 1 + - 2(3 cos9 with (3 = ^1- iM'^/s. The coefficients N, Ai and Bi are 
listed in Table 2, the effective couphng parameters Qi being defined by 



91^ Yl li^o 

a=L,R 



(6) 



92 



a=L,R 



ReKa{s) 



(7) 



^For a generalization ineluding anomalous couplings of vector states, sec Ref. [14]. The couplings to 
gauge bosons used here correspond to minimal couplings with k = 1 and A = 0. 
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^3 = 4 E (-] 

a=L,R \ ^ J 



with 



(8) 



(9) 



and \a denoting the Yukawa couphngs specified in Table 1. In the above, e = V^tto, 
= = 0, and Mz and Tz are the mass and the width of the Z boson, respectively. 



The electron and leptoquark electroweak charges are given by 



(10) 



CjnSq. 



where Q is the electric charge in units of e, T3 the third component of the weak isospin, 
and (cyj) the sine (cosine) of the weak mixing angle. 
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Table 2: Coefficients appearing in the cross section formula Eq. (5). 



The total cross section can be written in the following form [13, 14]: 



ttq; 



a — 



N[F,'gi + F^g2 + Figs) 



(11) 
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using the same effective couplings gi, i — 1, 2, 3 defined in Eqs. (6, 7, 8) and 




1, 

3fl + (3' (1-/3^)^ 



i5 



4[ 2/53 I-P)' 




(12) 



1 




1 + /3 



1-/3' 



The normalization factor is as given in Table 2. Asymptotically, for <^ s, the s- 
channel and the s/t- interference contributions scale like 1/s for scalar leptoquarks and ap- 
proach constant values for vectors. The pure ^-channel contribution scales with ^ ln(s/M^) 
for scalars, but grows like sjM^ for vectors. The latter behavior eventually leads to uni- 
tarity violation and thus indicates that the effective Lagrangian from which the above 
cross section formulae have been derived has to be embedded into a more fundamental 
theory at high energies. However, for Yukawa couplings obeying the experimental bounds 
of Table 1 this problem only arises far beyond the energy range considered here. 

At high-luminosity e^e~ colliders, beamstrahlung leads to a significant energy spread 
in the electron and positron beams. This effect can be described with the help of a radiator 
function Ge+e-(2^+) '^-1 *) where x± is related to the fractional energy loss, E± — x±^/s/2. 
The transverse momentum spread of the beam due to beamstrahlung can be neglected. 
Therefore, the cross sections including the effect of beamstrahlung can be obtained from 
Eqs. (5) and (11) by a simple convolution: 



Because of singularities at the end-points a;+ = 1 and = 1, Ge+e-{x+,X-) is split 
into 5-function terms (the probability for no energy loss by beamstrahlung) and a smooth 
function of x+ and X-. The radiator function depends on machine parameters and is 
numerically provided by the program circe [15]. 

The most important QED correction arises from initial state radiation (ISR). The 
universal part of this effect can again be described by a radiator function^: 




(13) 



= - .) [1 + _ ,„ _ (2 in. + -)J + .(1 - . - .)-^ 1,> (14) 



^£>e/e is actually known up to 0{a^), but since the non-universal contributions to 0{a) are not known 
for the process under consideration, we also do not include the O(a^) effects in -Dg/e- 
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Figure 1: Distribution of the production cross section in the effective collision energy \/s 
for the scalar leptoquark ^^^^Si. The three histograms refer to the Tesla (dashed), S-Band 
(full) and X-Band (dotted) options for y/s = 500 GeV, M ^ 200 GeV, X^^Xr^ 0. 



where E is the beam energy. The parameter e is a photon energy cutoff separating soft 
from hard photons. The (5-function term includes soft bremsstrahlung where photons 
have an energy < eE. The second term proportional to the ^^-function describes the 
emission of photons with E^ > eE. Similarly as in Eq. (13), the corrected cross section 
follows by convolution: 

da = dzDe/e{z)da{zE+, EJ) + dzDe/e{z)da{E+, zE^). (15) 

In addition, for s-channel diagrams, higher-order corrections due to vacuum polarization 
can be included by using the running fine-structure constant [16]. 

If one wants to take into account beamstrahlung simultaneously with ISR, one has to 
substitute in Eq. (15) da by da^g given in Eq. (13) and £'+ (-E-) by a;+y^/2 
The reduced center-of-mass energy squared after beamstrahlung and ISR is given by s = 

zs. The separation of Gg+e- ^"^^ ^eje into 5-function terms and smooth functions and 

the distinction of initial-state radiation from electrons and positrons requires to compute 
eight separate contributions to the cross sections. Technically, the mapping z ln(l — z) 
leads to numerically stable integration. 

Presently, three different options for the beam dynamics of an e^e' collider are being 
discussed. The differences in the distributions of the cross sections in s are exemplified in 
Fig. 1 for the leptoquark ~'^I'^S\ having the largest cross section. For the present study, 
we choose the Tesla option. The bulk of the events is accompanied by soft radiation with 
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a maximum at 1 — s/s ~ 10~^, but with a tail extending to a few percent. This hmits 

the possibihty to reconstruct the kinematics of final states with missing momentum by 
imposing the constraints of energy-momentum conservation. 

QCD corrections are not taken into account in our calculations. They are known only 
for the total cross section of scalar leptoquark production [17]. The correction factor 
5s — 1 + asF{(3) is shown in Fig. 2. At (3 < 0.2 the corrections exceed 100 %. In this 
region close to the threshold it will certainly be necessary to take into account bound 
state dynamics. At large /3, one expects corrections of the order of 20 % which increase 
the cross section. 




0.2 0.4 0.6 0.8 1 

/3 

Figure 2: QCD correction factor ds for the total production cross section of scalar lepto- 
quarks as a function of 13. 

In Fig. 3, finally, we illustrate the effects due to beamstrahlung and ISR on the inte- 
grated cross sections for representative scalar and vector leptoquarks. Except for scalars 
at very high y^, the total cross section is reduced typically by 10%: 1 to 2 fb for scalars 
and 60 fb for vectors. The reduction is particularly important for the discovery reach near 
threshold. 

Figure 4 shows the total cross section for leptoquark pair-production at two center- 
of-mass energies and for those species of scalar and vector leptoquarks which have the 
minimal and maximal cross sections in each class. From the figure it is clear that, given 
the mass, the measurement of the total cross section already provides an important piece 
of information on the type of the leptoquark produced. As demonstrated later, another 
important observable is the angular distribution. 

Additional information can be inferred from decay properties. They arc determined 
by the Yukawa couplings. The partial decay widths per channel (see Table 1) are given 
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0.014 n 1 1 1 1 1 n 1 rr 




400 600 600 700 800 900 1000 400 600 600 700 800 900 1000 

(a) (b) 

Figure 3: Effects of heamstrahlung and initial state radiation on the pair-production cross 
sections for the ^""^^^^o (a) and "^^^Vq (b) leptoquarks. The dotted curves are the Born 
cross sections, the thin full curves include heamstrahlung, and the thick full curves give 
the total cross section including both heamstrahlung and ISR. 



by [1] 

= t^\%rM for scalars, 

= ^^^l,R^ vectors. 

Quantitatively, for M = 200 GeV and Xl,r = e one has Tg (Ty) = 116 (77) MeV, i.e., 
very narrow states. Furthermore, for the leptoquarks considered here, the branching 
ratios Beq for charged lepton channels and Bi,q for neutrino channels add up to unity, 
Beq + Byq = 1, with Beq as givcu in Table 3. If one assumes that the members of a 
given isomultiplet of leptoquarks are almost mass-degenerate, one sees that the strongest 
bounds, M ^ 225 GeV for scalars and M ^ 300 GeV for vectors, obtained at the Tevatron 
for Beq — 1 apply to all leptoquark species except the two singlets ~^^^So and ~^^^Vo. The 
latter two may have branching fractions 1/2 < Beq < 1 so that only the weaker limits 
> 176 GeV and Mv;, > 270 GeV hold. Note that for first generation leptoquarks in 
this mass range, B^q must be very close to 1/2 (r = 0) or 1 (r = oo) unless Xl^r ^ 1- 
In more general scenarios where leptoquarks have additional decay channels, as is the 
case for example for squarks in super symmetric models with i?-parity violation, Beq and 
B^q are not fixed by the Yukawa couplings alone. They rather can be considered as 
independent parameters leading to more model-dependent mass bounds. We also mention 
the possibility for generic leptoquarks to have tiny Yukawa couplings and hence being very 
long-lived. In this case, they would have to form bound states (leptohadrons) which do 
not decay inside the detector. This requires special search strategies, e.g., for fractional 
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MiGeV 

Figure 4: Total cross sections for leptoquark pair-production at fixed center- of-mass en- 
ergies as a function of the leptoquark mass M assuming vanishing Yukawa couplings and 
including corrections due to beamstrahlung and ISR. 

charged particles. In our calculations, the branching ratios are always taken for generic 
leptoquarks in accordance with the Yukawa couplings chosen in the production process. 
Of course, sufficiently small \l,r has negligible influence on the production rates. 

The formulae and prescriptions summarized in this section are implemented in the 
Monte Carlo event generator LQPAIR^. This program simulates production and decay of 
leptoquarks and is interfaced to parton shower and hadronization routines of JETSET. In 
the case of initial state radiation, events containing a bremsstrahlung photon are generated 
according to Eq. (15). The logarithm In (AE'^/ml) in Eq. (14) arises from the integration 
over transverse momenta of the emitted bremsstrahlung photon or, equivalently, over the 

■''The program was developed from an earlier generator LQ2 by D. Gingrich [18]. The code is available 
from http : //www . desy . de/~ hspiesb/lqpair . html. 
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Isomultiplets 


Beq 


-'/'So 


1 + r 

2 + r 


-'/'So 
-'/'Vo 


1 


+2/3^1, -1/3^^^ -4/3^^ 

+y^V,, -2/3^^1, -V3y^ 


0, i 1 


-y'S^/„ -'/^S,„ 

-'/'v,„, -'/'v,„ 


" ,1 
1 + r' 


+'/^S,/„ -'/'Sy, 
^'/'Vl/2, -'/'%2 


0, 1 



Table 3: Branching ratios for charged lepton channels of generic leptoquarks (r — A|j/A|j. 



emission angle 9^ of the photon with respect to the incoming electron (positron) beam: 

4£;2 .+1 E , , 

In — —= \ dcosO^— — (17) 

where E 3> me is assumed. Therefore, transverse momenta are generated according to 
the integrand in Eq. (17). 



3 Experimental Considerations 
3.1 Detector Models 

In addition to beam- and bremsstrahlung, we shall take into account detector effects such 
as limited acceptance and resolution. Geometrical acceptance cuts (beam hole) reduce 
the cross section, in particular in the forward and backward regions, which is important 
in cases where large i-channel contributions (due to large Yukawa couplings) are present. 
Furthermore, measurement uncertainties leading to a smearing of the lepton and hadron 
four-momenta deteriorate the reconstruction of jets and thus make the reconstruction of 
leptoquark masses less certain. 

For our study we use two simple detector models interfaced to the leptoquark event 
generator. The main components are a tracking device and electromagnetic and hadronic 
calorimetry. Detector effects are described by Gaussian smearing of particle four-momenta. 
The single-particle detection efficiency is set to 98% and the threshold energy to 0.3, 0.15, 
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and 0.6 GeV for electromagnetic particles, charged hadrons, and neutral hadrons, respec- 
tively. Table 4 lists the main detector parameters for the two models. A more detailed 
discussion of the detector models will be given elsewhere [19]. Our final results refer 
to the dedicated 1 TeV detector. In selected cases we also compare the latter with an 
LEP/SLC-type detector. 





1 TeV detector 


LEP/SLC detector 


Electromagnetic calorimeter 

rr( F\ A 


Am = 0.10, = 0.01 


= 0.20, B,^ = 0.01 


Hadronic calorimeter 
E - VE^^"-' 


Aad = 0.50, Ehad = 0.04 


Ahad = 0.90, Ehad = 0.02 


Tracking (angular dependent) 
2 — 

Px,y 


Pr(90°) = 2 • 10-^ 


Pr(90°) = 6 • 10-^ 


Multiple scattering 

(^{Px,y,z) _ PmS 

Pl,y,z Px,y,z\/ sine 


Pms = 0.0015 


Pms = 0.0050 


Beam hole 


8.1° 


15° 



Table 4: Detector properties. 



3.2 Search Strategies 

Leptoquarks decay into a lepton, either a charged one or a neutrino, and a quark which 
is observed as a jet of hadrons. Therefore the final states to be searched for are 

(I) £+r + 2jets, (18) 

(II) £± + 2jets + p^iss, (19) 

(III) 2jets + p^iss (20) 

where the lepton and quark flavors are assumed to belong to one given generation. We 
shall focus on the first generation, although most of our results should also hold for second 
generation leptoquarks. 

Search I is straightforward. One has to identify two charged leptons in the event. The 
events should contain enough hadronic energy in order to allow the use of a jet algorithm 
to group the final state hadrons into two well-defined jets. Next, one combines one of the 
leptons with one of the jets. That combination which gives the smallest difference in the 
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invariant lepton-jet masses is accepted and the lepton-jet masses are associated with the 
leptoquark mass. 

Search II requires the identification of events with one charged lepton, two jets, missing 
momentum and the determination of Pmiss- Assuming the observed missing momentum 
to be carried by one single neutrino, the event analysis can be performed in complete 
analogy to search I . 

In search III, the signal is characterized by missing momentum and missing mass 
carried away by the unobserved neutrino pair. This is the most difficult case. Obviously, 
the measurements which can be carried out are not sufficient to reconstruct the leptoquark 
mass. 

For clean event identification and precise mass reconstruction we require the following 
cuts on energies and transverse momenta in search I to III: 

(I) two charged leptons with transverse momentum pip > 20 GeV, missing transverse 
momentum p^^^^ < 25 GeV, two jets with energies Ej > 10 GeV, and total visible 
energy E^ig > 0.9^/s■, 

(II) one, and only one, charged lepton with pip > 20 GeV, pjp'^^^ > 25 GeV, two jets with 
Ej > 10 GeV, and E^i^ > 0.6^^; 

(III) no charged lepton with pp > 20 GeV, pp'^^^ > 25 GeV, two jets with transverse 
momentum pp > 75 GeV each, and total hadronic energy Ehad < 300 GeV. The 
upper limit on Ehad corresponds to a lower limit for the missing mass of the neutrino 
pair. 

These general requirements are supplemented by additional cuts to suppress the back- 
ground processes which are discussed in the next subsection. 



3.3 Background Processes 

Vector boson pair-production is an obvious source of final states of the kind Eqs. (18) to 
(20) within the standard model: 

e+e- ^ ZZ, (21) 

e+e- ^ W+W- (22) 

with one Z (IF) decaying into a lepton pair and the other one into two jets. More 
generally, one has to consider four-fermion production 

e+e" hhqmi (23) 

comprising also single-resonant processes like e'^e" — > Zvu, ZU, Wiu with Z ^ qq and 
W — > qiq2, as well as completely non-resonant processes. Whereas on-shell W and Z 
production can be suppressed efficiently by rejecting events with pairs of final state leptons 
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or jets having a mass close to M\y or Mz, suppression of off-shell production requires 
additional cuts. We implement them as cuts on the various invariant masses Mi.ji^. 

Apart from the irreducible background due to final states containing the same particles 
as the signal events, there is in addition some reducible background due to processes with 
the same visible particles produced by different final states. In particular, 6-quarks and 
T-leptons have some probability to produce final states with an energetic electron, that 
may mimic the decay products of a leptoquark. 

A second important source of background is tt production: 

e+e~ tt (24) 
with the top-quarks decaying into b and W and the latter decaying leptonically. 

The cuts devised to suppress these background processes are summarized below. The 
indices ji and j2 label the two jets ordered by their energies, i.e., Ej^ > Ej^. h and I2 
denote the leptons ordered in such a way that 

\Mi,j, - Mi,j,\ < \Mi^j, - Mi,j,\ . (25) 

With this convention, the pairs liji and Z2j2 are candidates for the decay products of a 
leptoquark. In search (I) to (III) we require: 

(I) \Mij^ -Mz\> 10 GeV and \Mj^j^ - Mw,z\ > 10 GeV; 
Mi^i^ > 20 GeV, Mi^j^ > 20 GeV, Mi^j^ > 20 GeV. 

(11) \Mi^i^ -Mw\> 20 GeV and \Mj^j^ - Mw,z\ > 20 GeV; 
Mi^i^ > 20 GeV, Mi^j^ > 20 GeV, Mi^j^ > 20 GeV; 
Ehad > 150 GeV. 

(Ill) \Mi^i^ -Mz\>20 GeV and \Mj^j^ - Mw,z\ > 20 GeV; 
Mj^j^ < 400 GeV. 

Table 5 indicates the number of remaining events from various reducible and irreducible 
background processes. The estimates for four-fermion final states have been obtained with 
the help of WPHACT [20], those for tt production with PYTHIA [21]. Beamstrahlung and 
ISR are not included here. These corrections should be negligible for the present purpose. 

In search I, the dominant background is due to on- and off-shell boson pair-production 
with in the final state, i.e., e+e" —>■ e^e~qq and e+e" e^uqq'. The rate is much 
higher than at LEP2 energies estimated in a previous study of leptoquark production [22]. 
Whereas tt production is not a problem for search I, it makes the search for leptoquarks 
in channel II very difficult. In search III, the number of background events surviving the 
cuts is very large, since the final state cannot be reconstructed completely. In particular, 
final states with 6-quarks and r's contribute a large fraction. It seems unlikely that in 
this channel one can expect more than a consistency check of the searches I and II. 

Due to their huge cross sections, two-photon processes have also to be considered as a 
potential source of background although four-fermion final states only emerge at higher 
orders. We checked for a sample of 2 • 10^ 77 events generated with the help of the event 
generator PYTHIA (version 5.722) that no event passed the cuts for searches I and II, 
whereas in search III a negligibly small number of events remained. 
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e~VeCS + c.c. 


1 1 Q 

11. y 


141.0 


zoo 


1 ^ n 

iO.U 


oy .o 


iOOU 




— 


1.4 


495 


— 


1.3 


2410 


T qq 


0.4 


1.3 


888 


0.6 


5.1 


320 


other 4i 




0.1 


371 




8.0 


1250 


sum of 4f 


12.3 


17.6 


2010 


15.6 


53.7 


5540 




0.7 


132 


527 


1.5 


91.6 


97.8 


total 


13.0 


150 


2537 


17.1 


145.3 


5630 



Table 5: Number of background events surviving the cuts in the leptoquark searches (I) to 
(III) at^^ 500 GeV (C^ 20 fb'^) and 800 GeV (C^ 50 fb'^). 



4 Results 

In the following, we restrict ourselves to generic leptoquarks of the first generation which 
decay into or i/g and a jet. Except where stated differently we assume 

\l < 1, Ar < 1, (26) 

so that the Yukawa couplings have negligible effects on the production cross sections. 
Note that in this case the production rates are actually independent of the generation 
quantum number. However, the decay products of third generation leptoquarks contain 
the heavy flavors r, b and t which in turn decay into lighter leptons and jets. Hence, the 
search strategies explained in section 3 are not appropriate. 



4.1 Signal Distributions 

Figure 5 illustrates some characteristic distributions generated from the class II process 
e+e- ^ -^/^Si + ^/^Si e± ^ jj for M = 250 GeV and ^/s = 800 GeV in a high 
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Figure 5: Event distributions for the process e~^e~ — > ^^/"^Si + ^^'^Si — >• jj for 

M — 250 GeV at y/s — 800 GeV (search II). Masses are given in GeV and cross sections 
in ph. The definition of kinematic variables is explained in subsection 3.3. 
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statistics Monte Carlo run. Beamstrahlung, ISR, detector effects and cuts are taken 
into account as described above. In Fig. 5a one sees tlic sin^ ^-distribution cliaractcristic 
for the production of scalar particles which can nicely be reconstructed under the given 
assumptions. Figure 5b displays the distribution in the transverse momentum of the decay 
electron (positron), which is practically the same for the missing transverse momentum 
and the jet transverse momentum. The invariant-mass distributions of the lepton and 
the jet pairs shown in Figs. 5c and 5d have their maximum at about 300 GeV, while the 
leptoquark signal is clearly visible in the distribution of Mi^j^ of Fig. 5e. The tails in 
this distribution extending to smaller and higher masses are small. From Fig. 5f we see 
that the identification of the correct combination of leptons and jets is not perfect: the 
mass distribution of the 'wrong' combination shows some excess of events in the region of 
Mi^j^ ~ 250 GeV = M. 

4.2 Detection Efficiencies 

The cuts described in subsections 3.2 and 3.3 are considered a reasonable compromise 
between the two conflicting requirements of background suppression and preservation of 
the signal rate. It is certainly possible to further optimize them. Particularly worthwhile 
may be some fine-tuning for different ranges of center-of-mass energies and leptoquark 
masses. 

Table 6 gives an overview of the potential in the leptoquark searches I and II. Col- 
umn 2 displays the number of events observed divided by the respective branching ratios 
Bj — and Bu — 2BeqB^g or, equivalently, the number of leptoquarks produced and 
surviving the cuts in channel I and II, respectively. One can see that in both cases the 
signal efficiency is reasonably large, between 20 % and 40 %. There is no significant differ- 
ence in the efficiencies for scalar and vector leptoquarks. Furthermore, we find that for a 
LEP/SLD-type detector the reconstruction efficiency is lower by typically 20 % relative to 
the dedicated 1 TeV detector, depending on the channel and the leptoquark mass. This 
reduction originates mainly from the different size of the beam- holes (see Table 4) . 

The results in Table 6 are obtained for neghgibly small Yukawa couplings, i.e., for 
leptoquarks with very small natural widths (see Eq. (16)). Beamstrahlung, initial state 
radiation, and in particular hadronization effects lead to a considerable broadening of 
the observed leptoquark mass distribution. In order to estimate the mass resolution, we 
fitted the reconstructed mass distribution to a Gaussian. Column 4 of the table shows 
the standard deviation of such a fit. Note that the precision with which the peak position 
of the reconstructed mass distribution can be determined is, in most cases, better by a 
factor of 5 to 10, i.e., 6M = 0.5 to 0.8 GeV. From column 4 and Eq. (16) one clearly 
sees that a measurement of the mass distribution of leptoquark decay products will not 
allow to determine the intrinsic leptoquark width for accessible masses and values of A^^r 
allowed by the existing bounds. 

Finally, Table 6 reveals a systematic shift of the reconstructed leptoquark mass. Again, 
this is dominantly due to hadronization effects and depends strongly on the algorithm used 
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I 

M = 225 GeV 
^ = 500 GeV 


number of events 

-^events/ Bj 

{£ = 20 fb"^) 


signal efficiency (%) 


standard 
deviation 
(GeV) 


mass shift 
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II 

M = 225 GeV 
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^events/ Bji 

{£ = 20 fb-1) 








-'/'So 
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24.0 / 19.5 


6.5 / 8.0 


+1.3 / +1.6 


-''/'Vo 


269 / 234 


23.4 / 20.4 


6.7 / 7.3 


+1.6 / +0.9 


II 

M = 350 GeV 
= 800 GeV 


^events/ Bjj 

(£ = 50 fb-^) 








-'/'So 


16 / 12 


29.7 / 23.3 


5.9 / 6.7 


-0.8 / -1.3 


-''/'Vo 


488 / 397 


28.2 / 23.0 


5.8 / 6.9 


-0.7 / -1.5 



Table 6: Event number, signal efficiency, difference in reconstructed and nominal lep- 
toquark mass, and width of the reconstructed mass peak for two scenarios and in two 
search channels. The expected numbers of events in column 2 are divided by the respective 
branching ratios Bj = B^^, Bjj = 2BeqB^g. The first number given in each column refers 
to a dedicated 1 TeV detector, the second one to a LEP/SLC-type detector. 



for jet reconstruction (see for example [23]). We compared the JADE and Durham schemes 
with E, EO, P, and PO recombination and found that the Durham-P scheme leads to the 
smallest mass shifts. Other jet algorithms induce considerably larger shifts, for example, 
the JADE-E scheme a shift up to 10 GeV. Therefore, the Durham-P algorithm was used 
in the subsequent analysis. 

In channel 111, no mass reconstruction is possible. The signal efficiency ranges from 
40 to 45 % at ^ = 500 GeV and from 7 to 15 % at ^ = 800 GeV. 
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4.3 Sensitivity Limits 



From the above it becomes clear that the discovery potential in the three channels is 
very different. In other words, for a given leptoquark species, the observability strongly 
depends on the branching ratio for the decay into a charged lepton, here e^. As can be 
seen from Table 3, this branching ratio is fixed for all states with the exception of "^^^Sq, 
~^/^Vo, ~'^^^Si/2, and ~^^^Vi/2- In the latter cases, the branching ratio is determined by 
the ratio r = A^/A|. Below we shall present results for r — 0,1 and oo, corresponding to 
B^g = 1/2, 2/3, and 1. 

A first estimate of the sensitivity limits can be based on the total number of events. 
Requiring a 5(T effect, we determine the values of M for which the number of signal events 
is equal to or larger than b^N^g where N^g is the total number of background events. 
This is a sensible discovery criterion for search I where the number of background events 
is small. The results are collected in Table 7 for all leptoquark species of Table 1. In 
search I, scalar leptoquarks can be discovered for masses between 80 and 97% of \/s/2 
except for the states ^^^^Sq and ^^^^Si where only 70 % of ^/i/2 can be reached for purely 
left-handed couplings. For vector leptoquarks the mass reach is always larger than 95 % 
of ^/s/2. In searches II and III, the discovery limits are worse. In these channels, some 
scalars cannot be detected at all for masses above 100 GeV. It should also be noted that 
small, but non-negligible Yukawa couplings [13] and anomalous couplings [14] can lower 
the production cross sections and discovery limits in comparison to Table 7. 

A more refined analysis is possible if the number of events is sufficiently large to 
investigate the mass distributions. The signal is expected to stick out as a prominent 
peak over a flat background distribution. From an investigation of the mass distribution 
one would obtain a higher significance for leptoquarks with the masses shown in Table 
7 than the global 5a assumed there. This is demonstrated in Fig. 6a for a favorable 
case, while Fig. 6b illustrates a more difficult situation. In the latter, the 9 signal events 
expected for M = 200 GeV disappear in the ti background, whereas the 33 signal events 
for M = 160 GeV might still allow to identify an enhancement in the mass distribution. 
This, however, requires a mass resolution better than 5 GeV. 

Note that for Fig. 6 we have assumed = A/j <^ 1. Vanishing A^ forbids class II 
final states. For ~^/^5'i/2 (Fig. 6a) this is also the case for A^ = 0, while for "^^^^o the 
signal distributions shown in Fig. 6b have to be rescaled by the factor 8/9. For small, but 
non-negligible Yukawa couplings, the production rates may be either enhanced or reduced 
relative to the case of negligible Xl,r, whereas for sufficiently large couplings the t-channel 
contribution always leads to larger cross sections^. 

4.4 Leptoquark Couplings 

Assuming that the masses are known and the Yukawa couplings not too large, different 
species of leptoquarks can be distinguished ah^eady by their total production cross section 

**The influence of anomalous couplings is discussed in Ref. [14]. 
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Table 7: Discovery limits for leptoquarks (masses in GeV) at ^/s — 500 GeV (C — 20 
fb~^) and y/s — 800 GeV (C — 50 fb~^) requiring a 5a effect. Dashes indicate cases where 
the corresponding search is not possible, * means no sensitivity to masses above 100 GeV 
with the cuts considered. on 
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Figure 6: Signal and background distributions in the invariant mass Mi^j^ for the channel 
e^u+2jets and ^/s = 500 GeV, C ^ 20 fb''^: (a) ~^/^Si/2 production, (b) '^^^Sq production 
and M = 160 GeV (dashed), M = 180 GeV (dash-dotted) and M = 200 GeV (dotted). 
The full histograms show the dominant background from tt production. 



(see Fig. 4 and Table 1). The angular distribution gives an additional handle on the 
spin and the relative size of the couplings to gauge bosons and fermions. For most of the 
leptoquarks in Table 1 the s/i-channel interference is destructive reducing the cross section 
in the central region. By contrast, the pure t-channel contribution increases the cross 
section in the forward region. However, the Yukawa couplings have to be relatively large 
to make this an observable effect. In Fig. 7 the differential cross sections^ da/dO including 
all cuts, beamstrahlung and ISR are shown for ~'^/^Si and "^^^Vq which have comparable 
production cross sections (see Table 1). From the angular distribution it should not be 
difficult to distinguish scalars from vectors, but the shape of the distributions does not 
reflect the presence and size of Yukawa couplings as clearly as one would wish. 

Another difficulty follows from the fact that, in general, the detection efficiencies 
have an angular dependence which distorts the observable distributions, in particular in 
the forward and backward regions. Figure 8 illustrates the angular dependence of the 
detection efficiencies. For scalars (Fig. 8a) the cross section vanishes in the forward and 
backward region, leading to large statistical fluctuations. Hence, detector effects show up 
mainly in the normalization of the cross sections. However, for vector leptoquarks (Fig. 
8b), the efficiency decreases by almost 30% at small and large angles. This demands a 
precise unfolding of the efficiency, before one can hope to probe the Yukawa couplings. 

In order to evaluate the prospects for determining the couplings of scalars more quanti- 
tatively, we have fitted the measured angular distribution to the differential cross section 

^Hcrc, charge identification of the decay leptons is essential in order to remove the sign ambiguity of 
cos^. Experimentally this is possible for searches I and II even for very energetic electrons/positrons. 
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Figure 7: Angular dependence of the production cross sections for (a) ~^^^Si and (b) 
~^/^Vo for various Yukawa couplings {gL,9R) in units of e (M — 300 GeV, ^/s — 800 
GeV, C — 50fb~^, cuts included). 
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Figure 8: Angular dependence of the detection efficiency for a scalar (a) and a vector (b) 
leptoquark in channel I. The full histogram refers to the 1 TeV detector, the dashed one to 
a LEP/SLD-type detector. The dotted histogram is obtained with the detector simulation 
ofRef [24]. 
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da /d COS 9 given in Eq. (5) with the normahzation as a free parameter. Consequently, 
one can only determine ratios of the effective coupling parameters gi, g2, and g^, defined 
in Eqs. (6-8). Taking M = 200 GcV, ^ = 500 GcV and £ = 20fb-i, it turns out 
that for \g2/gi\ > 0.5 and ga/gi > 0.15 the ratios g2/gi and gs/gi can be determined to 
10% accuracy. The above ranges correspond to Yukawa couphngs Xl,r between 0.2 and 
0.6 (note that the relation between the g^ and Xl^r depends on the leptoquark quantum 
numbers). This result does not change strongly with s and M. For Xl,r < 0.1, the range 
allowed by low-energy data (see Table 1), it seems to be difficult to probe the Yukawa 
couplings by investigating the angular distributions. 

5 Conclusions 

Linear e'^e" colliders provide unique tools to search for leptoquarks independently of the 
size of their Yukawa couplings to lepton-quark pairs. The relatively clean environment of 
e~^e~ annihilation allows to reconstruct leptoquark masses from the decay products also 
in the presence of beamstrahlung, QED initial state radiation, and hadronization effects. 
The discovery limits depend on the quantum numbers of the particular leptoquark and 
are estimated to reach M — (0.9 to 0.95)^/s/2 in the most favorable cases. 

The masses can be reconstructed with a precision of 0.5 to 0.8 GeV assuming the 
detector performance expected for a dedicated 1 TeV detector and not accounting for 
systematic mass shifts. Past experience at LEP has shown that for the analysis of real data 
better measuring accuracies can be achieved relative to early simulations. Measurements 
of the total cross sections in combination with an analysis of the angular distributions 
will allow to clearly distinguish scalar from vector leptoquarks. Yukawa couplings can 
be determined probably only if they are of a strength comparable to the electromagnetic 
coupling. 

In summary, searches for leptoquarks at e'^e~ linear colliders are complementary to 
pp coUisions which also probe the existence of these novel states independently of the 
unknown Yukawa couplings, but with less power in distinguishing states with different 
quantum numbers. On the other hand, Yukawa couplings can be probed by searching 
for virtual leptoquark effects in e+e~ — > qq and pp —>■ e'^e~ , and most efficiently in ep 
coUisions. 
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